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Abstract 

The  deformable  leading  edge  (DLE)  concept  to  im¬ 
prove  the  blade  capability  in  lift,  drag  and  pitching  mo¬ 
ments  has  been  investigated  for  the  purpose  of  meeting 
new  rotor  maneuverability  and  susceptibility  require¬ 
ments.  The  advantages  and  disadvantages  of  this  con¬ 
cept  have  been  carefully  examined  with  limited  compu¬ 
tational  and  experimental  results.  This  work  showed 
that  this  concept  achieves  a  substantial  improvement 
in  lift  capability  and  also  reduces  the  drag  and  pitch¬ 
ing  moment  at  the  same  time.  Effects  of  various  pa¬ 
rameters,  such  as  Reynolds  number,  reduced  frequency, 
mean  angle  of  oscillation,  and  airfoil  shape,  on  the  per¬ 
formance  of  these  airfoils  were  also  investigated. 

Introduction 

In  forward  flight,  helicopter  rotor  blades  experience 
an  angle  of  attack  change  from  a  relatively  small  value 
on  the  advancing  sides  to  a  relatively  large  value  on  the 
retreating  side.  Under  certain  flight  conditions,  the  ro¬ 
tor  encounters  dynamic  stall.  Dynamic  stall  limits  the 
performance  of  helicopters  by  causing  severe  vibratory 
loads  due  to  abrupt  changes  in  lift,  drag  and  pitching 
moment.  This  load  condition  has  been  avoided  in  the 
past  by  defining  a  flight  envelope  which  keeps  the  air¬ 
craft  from  operating  in  the  region  of  dynamic  stall. 


However,  helicopters  are  now  being  considered  for  tasks 
which  require  vehicle  designers  to  consider  a  dramatic 
increase  in  the  performance  envelope.  Thus,  innovative 
ideas  must  be  incorporated  in  the  blade  design  to  delay 
or  suppress  dynamic  stall.  (See  Ref.l  for  more  details) 

In  order  to  meet  these  requirements,  computational 
and  experimental  works3  have  been  initiated  at  AFDD 
to  evaluate  the  effectiveness  of  various  concepts  such 
as  airfoil  slots/slats3’4,  deformable  leading-edge,  up¬ 
per  surface  blowing  (constant  and  cyclic),  flaperon,  and 
blade  flaps.  In  this  paper,  only  the  deformable  leading- 
edge  concept  will  be  discussed.  The  objective  of  DLE 
(Fig.  1)  is  to  delay  or  suppress  dynamic  stall  by  droop¬ 
ing  the  leading-edge  at  a  high  angle-of-attack  and  mak¬ 
ing  the  flow  pass  smoothly  around  the  leading-edge. 
To  investigate  the  static  and  dynamic  characteristics 
of  DLE,  an  incompressible  Navier-Stokes  code  (named 
ZETA),  based  on  a  velocity-vorticity  formulation,  was 
modified  to  treat  a  deforming  airfoil.  For  the  experi¬ 
mental  part,  the  Army  water  tunnel  at  Ames  was  uti¬ 
lised  to  measure  the  quasi-steady  and  unsteady  loads 
on  the  airfoil  during  the  pitch  oscillations.  The  influ¬ 
ence  of  Reynolds  number,  airfoil  shape,  mean  angle  of 
oscillation,  droop  angle  and  center  of  rotation,  and  re¬ 
duced  frequency  on  the  DLE  performance  was  also  in¬ 
vestigated. 

Description  of  the  Experiment 
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This  study  was  conducted  in  the  4000-liter,  closed- 
circuit  water  tunnel  facility  at  the  Aeroflightdynamics 
Directorate,  Ames  Research  Center  (Fig.  3).  This  was 
a  particularly  suitable  facility  for  this  investigation  be¬ 
cause  of  the  ease  of  obtaining  definitive  visualisations 
of  the  flow  and  the  advantage  of  examining  on-line  the 
resultant  loads  on  the  airfoil.  The  technique  for  visual¬ 
ising  the  flow  was  based  on  the  release  of  a  fluorescing 
dye  through  an  orifice  located  at  the  leading  edge  of 
the  airfoil.  Loads  were  measured  directly  by  an  exter- 
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nal  apparatus  that  served  as  both  support  and  balance 
for  the  airfoil. 

Two  airfoils  were  selected  for  this  study.  The  first 
was  a  VR-12  having  a  two-dimensional  planform  of  10 
cm  (chord)  by  21  cm  (span).  The  second  was  a  modified 
VR-12  with  a  droop  angle  of  13°  (case  2513).  Both 
airfoils  are  shown  in  Fig.  3.  The  test  section  measures 
31  cm  (height)  by  21  cm  (width),  and  the  airfoil  was 
positioned  so  that  it  spanned  the  width  of  the  section  to 
within  0.015  cm  on  either  side.  The  main  body  of  the 
airfoil  was  cast  from  urethane  and  the  spar  was  made 
from  stainless  steel. 

The  spar  of  the  airfoil  extended  through  the  test- 
section  windows  and  was  supported  by  lift  and  drag 
transducers  on  both  sides  (Fig.  4).  One  end  of  the  spar 
was  adjoined  to  an  instrumented  drive  shaft  through  a 
torsionally  stiff  coupling  so  that  airfoil  incidence  could 
be  set  and  the  pitching  moment  measured.  Static  fric¬ 
tional  moments  imparted  by  the  support  bearings  and 
seals  were  also  measured  and  treated  as  load  tares. 
Only  quantities  relating  to  the  airfoil  were  electrically 
instrumented,  and  these  were  incidence,  lift  (both  sides), 
drag  (both  sides),  total  pitching  moment,  and  the  bear¬ 
ing  and  seal  moments  (both  sides).  After  amplification, 
the  signals  were  were  transmitted  to  a  data  acquisition 
system  where  they  were  digitised,  averaged,  and  stored 
for  later  processing.  It  is  estimated  that  airfoil  inci¬ 
dence  was  measured  to  an  accuracy  of  0.2®  during  the 
test.  Lift  and  drag  measurements  are  considered  to  be 
accurate  to  0.01  newtons  and  the  pitching  moments  ac¬ 
curate  to  0.002  newton-m. 

The  dye  was  illuminated  by  an  xenon  strobe  during 
single-frame  exposures  and  by  an  argon-ion  laser  dur¬ 
ing  video  exposures.  In  both  cases,  the  light  sheet  was 
about  2  cm  wide  and  was  directed  through  the  upper 
test-section  window  and  across  the  flow  above  the  air¬ 
foil.  The  dye  that  was  released  from  the  leading  edge 
of  the  airfoil  was  transported  downstream  by  the  fluid 
in  the  boundary  layer  and  wake,  thereby  enabling  the 
thickness  and  eventual  separation  of  the  boundary  layer 
to  be  observed. 

The  tunnel  was  operated  at  an  average  dynamic 
pressure  of  0.019  kg /cm1,  which  for  an  ambient  tem¬ 
perature  of  21®C  and  an  airfoil  chord  of  10  cm  gives  a 
Reynolds  number  of  200,000.  The  unsteady  measure¬ 
ments  of  lift,  drag,  pitching  moment  were  made  with 
a  —  15®  +  10®  sin(2*/t)  and  at  a  reduced  frequency  k. 
The  reduced  frequency,  k,  is  related  to  the  oscillating 
frequency,  /,  by  the  expression  k  =  w/c/**,,  where  c  is 
the  chord  length  and  is  the  free  stream  velocity. 


Description  of  the  Theory 


The  ZETA*  code  uses  the  vorticity  vector  as  the 
primary  variable.  By  doing  so,  the  computation  of  an 
unsteady  flow  field  can  be  partitioned  into  its  kinematic 
and  kinetic  parts.  In  the  kinematic  part,  the  velocity 
field  is  related  to  the  vorticity  field  at  any  instant  of 
time  by  the  continuity  equation  and  the  curl  of  the 
velocity.  They  are 
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II 

> 

(1) 

V  x  v  =  w 

(2) 

In  the  kinetic  part,  the  vorticity  is  computed  at  each 
time  step  by  solving  the  vorticity  transport  equation, 
which  is 


— —  =  —(v  ■  V)w  +  V  •  (i/Vw)  (3) 

or 

where  v  and  w  are  the  velocity  and  vorticity  vectors 
respectively;  r  denotes  the  time  and  v  represents  the 
kinematic  viscosity  of  the  fluid.  When  the  flow  is  tur¬ 
bulent  v  includes  the  eddy  viscosity,  and  in  this  case 
the  Baldwin-Lomax  model  (Ref.  0)  is  used.  Eqs.  (1) 
and  (2)  can  be  recast  into  an  integral  form  for  the  veloc¬ 
ity  vector  using  the  fundamental  solution  of  an  elliptic 
equation: 


J(f,  t)  =  -  /  «. 

Jr 


x  QdR„ 


+  J  [(®o  •  "<>)  -  (®„  x  ri0)x]QdB0  +  i 


(4) 


where  r  is  a  position  vector  and  is  the  free  stream 
velocity;  R  is  the  fluid  domain  contained  between  the 
far  stream  boundary  and  the  solid  surface  S;  n  is  a 
unit  normal  vector  on  5  directed  outward  into  R;  the 
subscript  “o”  indicates  that  a  variable  or  an  integra¬ 
tion  is  in  the  f„  space;  and  Q  is  the  gradient  of  the 
fundamental  solution.  For  two  dimensional  flows, 


Q  =  -  ® ~  r- 

2 *!?,  -  f]1 

It  should  be  noticed  from  Eqs.  (3)  and  (4)  that  the  com¬ 
putation  of  any  unsteady  flow  problem  can  be  confined 
to  the  vortical  region  (without  regard  for  the  remain¬ 
der  of  the  flow  field).  This  ability  offers  a  significant 
savings  in  both  computer  time  and  storage.  For  the 
present  study  of  dynamic  stall,  the  vortical  flow  region 
consists  of  the  attached,  the  separated  and  wake  sones. 
Whenever  the  flow  is  attached,  a  boundary  layer  solver 
is  employed  to  increase  the  computing  efficiency.  The 
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vorticity  values  on  the  solid  surface  require  special  con¬ 
sideration  since  the  non-slip  boundary  condition  on  the 
solid  surface  is  required  for  the  velocity  field.  As  shown 
in  Ref.  7,  the  boundary  vorticity  values  are  obtained 
evaluating  Eq.  (4)  on  S, 

»(»\s,r)  =  -  /  x  QdR0 
Js+ 

-  J  uax  QdR0  (6) 
+  J  K?o  •  «o)  -  (»„  X  HoJxJQdB,,  +  V", 

where  F5  is  a  position  vector  on  the  solid  surface  5 
where  the  vorticity  is  unknown.  The  domain  integral 
over  R  in  Eq.  (4)  is  divided  into  two  parts:  one  is 
the  integral  over  an  infinitesimal  thin  layer  S+  adja¬ 
cent  to  5  containing  the  unknown  vorticity,  and  the 
other  is  the  integral  over  the  remaining  vortical  region 
R~ .  The  unknown  vorticity  within  the  infinitesimally 
thin  layer  is  represented  by  the  vortex  sheet  strength 
7.  The  strength  7  can  be  determined  from  Eq.  (6) 
if  the  vorticity  distribution  in  R~  and  velocity  on  S 
are  known.  The  conservation  of  the  total  vorticity  is 
needed  to  determine  the  strength  uniquely. 

The  numerical  procedure  incorporates  a  conformal 
mapping  to  transfer  the  physical  domain  onto  a  simple 
region  with  suitable  computation  grids.  The  initial  so¬ 
lution  starts  bom  a  potential  flow  solution  around  the 
solids  undergoing  an  impulsive  start  motion.  A  time 
marching  procedure  computes  the  initial  vortex  sheet 
strength  7  in  the  infinitesimally  thin  layer  surrounding 
the  solids.  The  velocity  field  corresponding  to  this  ini¬ 
tial  vorticity  field  is  computed  by  Eq.  (4).  The  vorticity 
is  then  diffused  from  the  thin  layer  and  convected  away 
by  the  velocity  field.  The  process  is  repeated  at  each 
time  step  until  either  a  converged  or  a  periodic  solution 
is  reached.  A  detailed  description  of  the  mathemati¬ 
cal  formulation,  conformal  transformation  and  solution 
procedure  is  given  in  Ref.  5.  The  procedure  to  solve 
the  variable-droop  leading-edge  airfoil  requires  the  con¬ 
formal  transformation  of  the  airfoil  coordinates  at  each 
time  step. 

Results  and  Discussion 


Effects  of  droop  angle  and  center  of  rotation 

The  deformable  leading-edge  (DLE)  approach  to  de¬ 
laying  or  suppressing  the  dynamic  stall  is  based  on 


drooping  the  leading-edge  part  of  the  airfoil  at  high  an¬ 
gles  of  attack  so  that  the  flow  passes  smoothly  around 
the  leading-edge.  There  are  two  geometric  parameters 
describing  the  drooped  leading-edge;  one  is  the  center 
of  rotation(zc)  and  the  other  is  the  droop  angle(tf),  as 
illustrated  in  Fig.  1.  In  our  notation,  for  example,  case 
2513  means  that  the  center  of  rotation  is  located  at  25 
%  of  the  chord  and  the  droop  angle  is  13s.  A  series  of 
numerical  tests  of  the  static  and  dynamic  performances 
of  drooped  VR-12  airfoils  was  carried  out,  with  the  cen¬ 
ter  of  rotation  and  the  droop  angle  as  parameters,  at  a 
Reynolds  number  of  200,000  and  a  reduced  frequency 
of  0.1.  Dynamic  lift,  drag,  and  pitching  moments  dur¬ 
ing  an  oscillation,  with  the  droop  angle  and  the  center 
of  rotation  as  parameters,  are  shown  in  Fig.  5  and  Fig. 
6.  The  results  indicate  the  following: 

(1)  As  the  droop  angle  increases  with  a  fixed  cen¬ 
ter  of  rotation,  the  stall  angle  and  the  maximum  lift 
decreases,  but  the  drag  and  pitching  moments  are  also 
greatly  reduced.  When  the  droop  angle  increases  be¬ 
yond  15s,  the  drag  and  pitching  moment  are  not  no 
longer  reduced  and  the  stall  angle  continues  to  decrease. 
(See  Fig.5) 

(2)  As  the  center  of  rotation  moves  rearward  (ze  = 
.15  — ►  .25)  with  a  fixed  droop  angle  (6  =  15°),  the 
stall  angle  increases  and  maximum  drag  and  pitching 
moments  are  reduced.  (See  Fig.  6) 

(3)  Based  on  a  variation  of  the  droop  angle  and  the 
center  of  rotation,  case  2513  appears  to  be  superior 
to  the  no  droop  case,  offering  a  reduction  in  drag  and 
pitching  moment  and  an  increase  in  lift.  (See  Fig.  7) 
Note  that  case  2513  also  turned  out  to  be  the  best  case 
for  Re=2,000,000  following  a  similar  parametric  study, 
(see  Lee  and  McAlister*  for  details) 

Water  tunnel  measurements  of  the  basic  VR-12  and 
a  drooped  VR-12  (case  2513)  airfoil  are  compared  with 
the  computed  results  in  Fig.  7.  La  the  experiment, 
case  2513  showed  a  significant  decrease  in  lift  hystere¬ 
sis  with  about  40  %  of  reduction  in  the  magnitude  of 
the  drag  and  pitching  moment.  The  gain  in  overall  lift 
(integrated  over  the  cycle  of  oscillation)  was  achieved 
mainly  on  the  downstroke  since  the  stall  is  delayed  and 
attenuated  by  drooping  the  leading  edge.  On  the  other 
hand,  the  computed  results  predicted  earlier  stalls  in 
both  cases  (especially  in  a  drooped  case),  and  conse¬ 
quently  overestimated  the  drag  and  pitching  moments 
reductions. 

Verification  of  computational  results 

In  this  section,  the  experimental  results  from  the 
water  tunnel  test  are  compared  with  the  predicted  re- 
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stilts.  Fig.  8  shows  a  comparison  of  experimental  and 
computational  flow  visualisations  of  the  basic  VR-12 
airfoil  at  conditions  of  Re  =  200, 000., a  =  15°+10°  sin  ait, 
and  k  =  0.1.  Up  to  15°  on  the  upstroke  the  flow  re¬ 
mains  attached  in  both  cases.  Starting  at  an  angle  of 
20° ,  the  experiment  shows  a  much  larger  separated  re¬ 
gion  than  does  the  computation. 

The  pitching  moment  hysteresis,  rather  than  the  lift 
hysteresis,  might  be  a  more  reasonable  measure  of  the 
prediction  capability  of  a  code0.  Following  this  point 
of  view,  the  code  seems  not  to  predict  the  stall  onset 
accurately,  but  the  overall  accuracy  is  acceptable  (Fig. 
7).  One  place  that  is  easy  to  blame  for  the  discrepancy 
is  turbulence  modeling.  Needless  to  say,  turbulence 
modeling  should  be  improved  to  simulate  the  dynamic 
stall  phenomenon  accurately.  Unfortunately  none  of 
the  eddy-viscosity  turbulence  models  (from  algebraic 
to  two-equation)  can  predict  separation  accurately  in 
all  cases10 .  Another  possible  reason  is  transition  mod¬ 
eling.  This  is  particularly  true  if  the  Reynolds  number 
is  below  500,000  and  a  separation  is  followed  by  a  reat¬ 
tachment  of  the  flow,  forming  what  is  known  as  a  “sepa¬ 
ration  bubble”.  This  phenomenon  can  not  be  properly 
simulated  without  an  accurate  transition  model.  As 
an  example,  Fig. 9  shows  a  comparison  between  experi¬ 
mental  and  computational  results  for  a  very  slow  oscil¬ 
lation  (quasi-steady).  Measurement  of  the  stall  onset 
reveals  a  significant  discrepancy  from  the  calculated  re¬ 
sults  wherein  the  flow  was  set  to  be  fully  turbulent.  In 
many  cases,  experimentalists  will  trip  the  airfoil  to  in¬ 
duce  artificial  transition  at  a  pre-determined  position, 
which  makes  it  easier  to  compare  with  computed  re¬ 
sults  when  the  code  does  not  have  a  transition  model. 

A  series  of  experiments  for  the  tripped  airfoils,  as  a  ex¬ 
tension  of  the  present  study,  is  in  progress  and  will  be 
shown  in  Ref.  8. 

Deformable  leading  edge  at  a  higher  Re 

Since  a  fixed-droop  airfoil(case  2513)  showed  a  sig¬ 
nificant  enhancement  in  lift,  drag,  and  pitching  moment 
hysteresis  at  a  low  Reynolds  number  (Re  =  200, 000), 
it  is  worth  while  to  investigate  the  performance  at  a 
higher  Reynolds  number  (Re  =  2,000,000).  For  this 
investigation  at  the  present  stage,  the  ZETA  code  was 
used  to  explore  two  situations.  One  is  a  fixed-droop 
leading-edge  airfoil  (case  2513)  and  the  other  is  a  variable- 
droop  leading-edge  airfoil  (YDLE  hereafter).  The  VDLE 
concept  is  desirable  for  keeping  the  original  shape  until 
the  airfoil  meets  an  angle  of  attack  slightly  below  the 
stall  angle  in  order  to  obtain  as  high  a  lift  as  possible 
cu  the  upstroke.  In  this  section,  the  detailed  results 


for  a  fixed-droop  (case  2513)  and  a  VDLE  airfoil  are 
described.  The  conditions  are  Re  =  2, 000, 000.,  o  = 
15s  +  10°  sin  cut,  and  k  =  0.1.  For  the  VDLE  case, 
the  droop  angle,  0,  versus  the  angle-of-attack,  a ,  was 
selected  as  follows: 

(  9  =  a  -  13.0  if  a>  13.0 

\  0  =  0.0  if  a  <  13.0 

Before  proceeding  with  the  detailed  calculation,  it 
should  be  verified  that  the  code  works  properly  at  a 
Reynolds  number  of  2,000,000  since  the  computed  re¬ 
sults  revealed  a  poor  prediction  of  quasi-steady  results 
at  Re  =  200, 000.  In  Fig.  10,  the  predicted  results 
for  steady  angles-of-attack  are  compared  with  available 
experimental  data.11’10  For  the  NACA  0015  and  VR-7 
airfoil  cases,  the  computed  results  predict  the  stall  on¬ 
set  and  maximum  lift  coefficients  quite  well.  Note  that 
both  experiments  were  carried  out  at  a  Mach  number  of 
0.1,  where  compressibility  effects  should  be  negligible. 

The  static  lift,  drag  and  pitching  moment  were  cal¬ 
culated  and  measured  at  discrete  angles  from  0°  to  30°. 
The  static  and  dynamic  load  characteristics  of  the  ba¬ 
sic  VR-12  airfoil,  the  drooped  VR-12  (case  2513),  and 
a  VDLE  case  are  compared  in  Fig.  11.  The  static  lift 
coefficient  for  the  no-droop  case  simply  increases  un¬ 
til  a  =  15°.  Slightly  above  a  =  18°,  the  airfoil  stalls, 
which  causes  Ci  to  drop  to  the  minimum  value.  For 
case  2513,  the  airfoil  also  stalls  at  a  =  13°,  but  the 
magnitude  is  much  smaller  since  the  stall  occurs  near 
the  trailing-edge  instead  of  the  leading-edge.  This  air¬ 
foil  eventually  stalls  from  the  leading-edge  (at  a  =  25° ) 
and  the  lift  drops  quickly.  The  variable-droop  leading- 
edge  (VDLE)  case  performs  much  better  than  the  pre¬ 
vious  two  cases  in  the  static  mode.  A  comparison  of  the 
drag  obtained  in  these  three  cases  most  clearly  shows 
that  the  drooped  airfoil  delays  the  static  stall.  The  drag 
coefficient  for  the  basic  VR-12  shows  an  abrupt  increase 
at  a  =  16°  due  to  stall.  In  case  2513,  the  drag  increases 
more  gradually  and  reaches  a  smaller  magnitude  than 
obtained  in  the  basic  case.  For  the  VDLE  case,  the  drag 
is  generally  less  them  that  of  case  2513.  The  pitching 
moments  of  both  drooped  leading-edge  airfoils  have  less 
variation  with  angle-of-attack  than  does  the  basic  case. 

To  better  appreciate  the  extent  of  the  unsteady  hys¬ 
teresis  in  the  lift  curves,  the  static  lift  curve  (dotted 
line)  for  the  basic  VR-12  airfoil  has  been  included  in 
Fig.  11.  For  the  basic  VR12  airfoil,  the  dynamic  stall 
effect  is  quite  evident.  The  lift  increases  almost  linearly 
up  to  a  =  23s,  then  increases  rapidly  due  to  the  forma¬ 
tion  of  leading-edge  stall  vortex,  and  finally  experiences 
a  sudden  drop  due  to  the  shedding  of  the  leading-edge 
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vortex  into  the  downstream  wake.  Note  that  the  an¬ 
gle  of  the  leading-edge  vortex  formation  is  delayed  by 
2°,  as  compared  with  the  case  at  Re=  200,000.  When 
the  vortex  is  detached  from  the  airfoil,  the  airfoil  ex¬ 
periences  large  fluctuations  in  lift,  drag  and  pitching 
moment,  as  seen  in  Fig.  11.  The  lift  hysteresis  loop 
for  the  drooped  VR-12  (case  2513)  shows  that  it  does 
not  encounter  deep  dynamic  stall,  consequently  both 
the  drag  and  the  moment  hysteresis  loops  also  do  not 
experience  any  abrupt  increase,  as  compared  with  the 
basic  VR-12  airfoil.  The  VDLE  case  reveals  an  addi¬ 
tional  gain  in  lift  on  the  upstroke  and  more  reduction  in 
drag  and  pitching  moment.  Note  that  the  VDLE  case 
generally  performs  better  than  case  2513  in  the  range 
of  low  angles-of-attack. 

Computed  surface  pressure  distributions  and  stream¬ 
lines  at  selected  angles-of-  attack  during  the  oscillation 
for  the  basic  and  the  drooped  VR-12  (case  2513)  airfoil 
are  shown  in  Figs.  12  and  13.  The  basic  VR-12  expe¬ 
riences  a  strong  leading-edge  stall  at  a  =  25°  and  has 
lost  most  of  the  suction  peak.  On  the  other  hand,  the 
drooped  VR-12  (case  2513)  experiences  only  trailing- 
edge  stall  even  at  a  =  25° ,  and  consequently  maintains 
its  suction  peak.  A  comparison  of  vorticity  contours 
in  Fig.  14  indicates  that  the  basic  VR-12  induces  a 
large  leading-edge  separation  at  a  =  25°  whereas  the 
drooped  airfoil  does  not. 

Effect  of  reduced  frequency 

The  contrast  in  performance  of  the  basic  and  drooped 
VR-12  (case  2513)  airfoil,  with  the  reduced  frequency 
as  a  parameter,  is  shown  in  Fig.  15.  The  stall  an¬ 
gle  of  the  basic  airfoil  is  more  sensitive  to  changes  in 
reduced  frequency  than  that  of  the  drooped  airfoil,  a 
trend  which  is  clearly  evident  in  all  (lift,  drag  and  pitch¬ 
ing  moment)  hysteresis  curves.  The  magnitudes  of  the 
lift  and  drag  coefficients  are  increased  with  frequency 
in  both  cases.  This  implies  that  the  performance  is  en¬ 
hanced  by  drooping  the  leading  edge,  regardless  of  the 
reduced  frequency.  Note  that  in  the  quasi-steady  case 
(k  =  0.0020),  the  drooped  airfoil  does  not  exhibit  any 
hysteresis  in  the  lift,  drag,  or  pitching  moment  curves. 

Effect  of  mean  angle 

Fig.  10  shows  the  effect  of  the  mean  angle  of  os¬ 
cillation  on  the  unsteady  loads  of  both  the  basic  and 
the  drooped  cases.  For  the  three  cases  studied  (am  = 
10°,  15°,  and  20°),  the  drooped  airfoil  has  better  per¬ 
formance  in  lift,  drag,  and  pitching  moment.  However, 
the  improvements  obtained  for  am  =  10°  and  am  =  20° 


are  not  as  striking  as  those  obtained  for  am  =  15°.  The 
reason  is  that  the  case  2513  was  determined  as  the  best 
configuration  for  am  =  15°,  but  not  for  am  =  10° 
or  am  =  20°,  which  implies  that  the  performance  en¬ 
hancement  obtained  by  drooping  the  leading-edge  de¬ 
pends  on  the  mean  angle  to  some  extent. 

Effect  of  airfoil  shape 

To  investigate  how  a  different  airfoil  shape  might 
affect  dynamic  stall  after  drooping  the  leading  edge, 
the  SSC-A07  airfoil18  was  selected.  In  Fig.  17,  dy¬ 
namic  lift,  drag,  and  pitching  moments  are  compared 
for  the  basic  and  drooped  SSC-A07  airfoils.  The  over¬ 
all  tendency  is  very  similar  to  the  VR-12  case,  which 
showed  that  drag  and  pitching  moment  are  greatly  re¬ 
duced,  with  some  gain  in  lift,  due  to  a  reduction  in  the 
magnitude  of  the  dynamic  stall. 

Conclusions 

The  present  study  considered  the  effect  of  a  de¬ 
formable  leading-edge  for  dynamic  stall  control.  With 
limited  experimental  and  computational  data,  the  ma¬ 
jor  concluding  remarks  are  as  follows: 

(1)  By  drooping  the  leading  edge  during  a  pitch  os¬ 
cillation,  the  lift  performance  can  be  significantly  en¬ 
hanced  and  the  drag  and  the  pitching  moment  greatly 
reduced  at  both  low  (Re  =  200, 000)  and  high  (Re  = 
2, 000, 000)  Reynolds  numbers.  The  gain  in  lift  occurs 
mostly  during  the  downstroke,  rather  than  on  the  up¬ 
stroke,  as  a  result  of  an  attenuation  of  the  magnitude 
of  the  dynamic  stall.  The  deformable  leading-edge  also 
appears  to  be  more  effective  than  the  slatted  airfoil  (see 
Ref.l)  for  reducing  drag  and  pitching  moment  fluctua¬ 
tions. 

(2)  The  benefits  from  a  drooped  leading  edge  are  not 
sensitive  qualitatively  to  the  reduced  frequency  and  the 
airfoil  shape,  but  might  be  sensitive  to  the  mean  angle. 

(3)  The  ZETA  code  predicted  the  dynamic  stall 
phenomena  quite  well  qualitatively,  but  was  not  exact 
quantitatively.  The  code  predicted  the  stall  onset  rea¬ 
sonably  in  cases  of  high  Reynolds  number,  but  did  not 
in  the  cases  of  low  Reynolds  number  with  a  low  reduced 
frequency.  The  discrepancies  are  partially  attributed  to 
the  lack  of  appropriate  transition  and  turbulence  mod¬ 
els. 

(4)  As  an  extension  of  this  study,  the  compressibility 
effect  on  the  deformable  leading  edge  should  be  inves¬ 
tigated. 
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[Pig.  1]  Schematic  drawing  of  drooped  leading  edge  and  [Pig.  2]  Water  tunnel  at  NASA- Ames  Research  Center, 
its  coordinates. 


BASIC  VR-12 


DROOPED  VR-12  (CASE  2513) 


[Fig.  3]  Airfoil  shapes  of  basic  and  drooped  (case  2513)  VR-12. 
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(Fig.  4]  Drawing  of  the  test  section  of  water  tunnel,  showing  a  balance  system  and  model  installation. 
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VR-07  airfoil 


NACA  0015  airfoil 


[Fig.  10]  Comparison  of  measured  and  computed  steady  lift  coefficients  of  NACA  0015  and  VR-7  airfoils  at  a  higher 
Reynolds  number  (  Re  =  2, 000, 000.) 
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Basic  VR-12,  Ra  =  2,000,000.,  a  =  15*  +  10*  sin  (dt,  k  =  0.1 


o  -2  .4  .6  .8  1.0  0  .2  .4  .6  .8  1.0 

Chord  Chord 


[Fig.  12]  Evolution  of  the  surface  pressure  distribution  and  streamlines  of  the  basic  VR-12  airfoil  during  oscillation 
at  conditions  of  Re  =  2, 000, 000.,  a  =  15°  +  10°s*nwt,4  =  0.1. 
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Drooped  VR-12  (case  2513),  R«  =  2,000,000.,  a  =  15*  +  10*  sin  tot,  k  =  0.1 


0  .2  A  .6  .8  1.0  0  .2  .4  .6  .8  1.0 

Chord  Chord 


[Fig.  13]  Evolution  of  the  surface  pressure  distribution  and  streamlines  of  the  drooped  (case2513)  VR-12  airfoil 
during  oscillation  at  conditions  of  JZe  =  2, 000, 000.,  o  —  15®  +  lO'iitwi,  1  -  0.1. 
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